use "data.dta",clear
xtset cowcode year

*** 表-2

logit bri_join_dum neighjoin400_bri_lag1  initiative_cnty years_diplomatic_lag1 IdealPointDistance_chn_lag1 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 dep_export_lag1 logviolences_lag1 i.year , vce(cluster iso)
estimates store m1
meglm bri_join_dum neighjoin400_bri_lag1  initiative_cnty years_diplomatic_lag1 IdealPointDistance_chn_lag1 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 dep_export_lag1 logviolences_lag1  ||iso:,  family(binomial) link(logit)
estimates store m2
meglm bri_join_dum neighjoin400_bri_lag1  initiative_cnty years_diplomatic_lag1 IdealPointDistance_chn_lag1 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 dep_export_lag1 logviolences_lag1  ||year:,  family(binomial) link(logit)
estimates store m3


esttab m1 m2 m3 using "table2.rtf",se parentheses ///
	nonumbers mtitles("1" "2" "3") ///
  star(* 0.10 ** 0.05 *** 0.01)  replace 


* 附表-2
heckprobit bri_join_dum neighjoin400_bri_lag1 years_diplomatic_lag1 IdealPointDistance_chn_lag1 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 dep_export_lag1 logviolences_lag1 , select(initiative_cnty = years_diplomatic_lag1 IdealPointDistance_chn_lag1 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 dep_export_lag1 logviolences_lag1)  
estimates store s2
**vce(cluster iso)

meglm bri_join_dum neighjoin400_bri_lag1  initiative_cnty years_diplomatic_lag1 IdealPointDistance_chn_lag1 IdealPointDistance_usa_lag1  arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 dep_export_lag1 logviolences_lag1  ||iso:,  family(binomial) link(logit)
estimates store s4


esttab s2 s4 using "Tex/Tables/table4.rtf",se parentheses ///
	nonumbers mtitles("1" "2" "3" "4") ///
  star(* 0.10 ** 0.05 *** 0.01)  replace 


** robustness check

* 附表-3

meglm bri_join_dum  neighjoin400_bri_lag1  initiative_cnty years_diplomatic_lag1 IdealPointDistance_chn_lag1 ///
 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 ///
 dep_export_lag1 logviolences_lag1  sco_full_lag1 asean_memb_lag1  ||iso:,  family(binomial) link(logit)
estimates store r1

meglm bri_join_dum  neighjoin400_bri_lag1  initiative_cnty years_diplomatic_lag1 IdealPointDistance_chn_lag1 ///
 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 ///
 dep_export_lag1 logviolences_lag1   bits_chn_lag1 pta_china_lag1 partner_level_lag1 ||iso:,  family(binomial) link(logit)
estimates store r2

meglm bri_join_dum  neighjoin400_bri_lag1  initiative_cnty years_diplomatic_lag1 IdealPointDistance_chn_lag1 ///
 arms_logexportfromCHN_lag1 arms_logimporttoCHN_lag1 v2x_polyarchy_lag1 ctygdppc_log_lag1   dep_import_lag1 ///
 dep_export_lag1 logviolences_lag1  minidist_log_lag1 ||iso:,  family(binomial) link(logit)
estimates store r3

esttab r1 r2 r3 using "Tex/Tables/robust.rtf",se parentheses ///
	nonumbers mtitles("1" "2" "3") ///
  star(* 0.10 ** 0.05 *** 0.01)  replace 

